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*sO ! Abstract 

^ ■ We suggest that gluon- induced dissociation and screening of the T(nS) states 

ft, . explain the suppression of the T(2S'+3S') states relative to the T(1S) ground 

state that has been observed by CMS in PbPb collisions at a/J/vjv = 2.76 TeV 
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at the CERN LHC The minimum-bias gluodissociation cross sections of the 
1S — 3S states are calculated using a screened Cornell potential and a thermal 



! gluon distribution. The 3S state dissolves due to screening before sizeable 

O 1 

gluodissociation occurs, but for the 2S and IS" states there is an interplay 
between screening, gluodissociation, and feed-down from the Xfe(2P) and 
Xfe(l-P) states. Based on a schematic approach, we find that the calculated 
suppression of the T(2S) and T(3S) states relative to T(1S) is consistent 
with the CMS result, but allows for additional suppression mechanisms. The 
T(1S) suppression through gluodissociation is, however, in good agreement 
with the CMS data. 
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The suppression of quarkonium states is one of the most promising probes 
for the properties of the quark-gluon plasma (QGP) that is generated in 
heavy-ion collisions at high relativistic energies. In the QGP the confining 
potential of heavy quarkonium states is screened due to the interaction of the 
heavy quark and the antiquark with medium partons and hence, charmonium 
and bottomium states successively melt Q at sufficiently high temperatures 
Tdiss beyond the critical value T c ~ 170 MeV. 

Charmonium suppression has been studied since 1986 in great detail 
both theoretically, and experimentally at energies reached at the CERN Su- 

S;r Proton Synchrotron SPS, BNL Relativistic Heavy-Ion Collider (RHIC) 
, [3, 0, 5], and CERN LHC (§[ 7]. The precise origin is still under investi- 
gation, in particular at LHC energies where regeneration due to statistical 
recombination of c and c in the quark-gluon plasma could be relevant, coun- 
teracting the J/ty dissociation in central collisions and contributing to the 
measured nearly flat suppression factor as function of centrality for > 

a. 

Bottomium suppression is expected to be a cleaner probe. The T(IS') 
ground state with invariant mass 9.46 GeV is strongly bound, the threshold 
to BB decay is at 1.098 GeV. Its lifetime of 1.22-10~ 20 s is about 1.7 times 
as large as the one of J/ty(lS) in elementary collisions. It melts as the last 
quarkonium in the QGP (depending on the potential) only at 4.10 T c [8j, 
whereas the 2S (10.02 GeV) and 35" (10.36 GeV) states melt at about 1.6 
and 1.2 T c , respectively. Even at LHC energies the number of bottom quarks 
in the QGP remains small such that statistical regeneration of the T states 
is unimportant. 
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T suppression in heavy-ion collisions has recently been observed for the 
first time both by the STAR experiment at RHIC [9(, and by the Compact 
Muon Solenoid (CMS) experiment at LHC [lOj. The latter includes an obser- 
vation of the enhanced suppression of the 2S + 35 relative to the IS ground 
state, whereas the IS suppression itself is considered in by CMS. 

This result is most likely not due to differences in the direct bottomium 
production mechanism in pp vs. PbPb collisions since nuclear modification 
of the parton distribution functions (shadowing) should affect all three states 



in a similar fashion 
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In this Letter we investigate the suppression of T(15), (2S), (3S) states 
at LHC energies due to screening and gluon-induced dissociation, including 
feed-down from the ^(lP) and x fc (2P) states. Whereas gluodissociation 
below T c is not possible due to confinement, it does occur above T c where 
the color-octet state of a free quark and antiquark can propagate in the 
medium. The process is relevant below the dissociation temperature T diss 
that is due to Debye screening, and its significance increases substantially 
with the rising gluon density at LHC energies. 

In the midrapidity range \y\ < 2.4 where the CMS measurement [10| has 
been performed, the temperature and hence, the thermal gluon density is 
high, and causes a rapid dissociation in particular of the 2S and 3S states, 
but also of the IS ground state. At larger rapidities up to the beam value of 
Ubeam = 7.99 and correspondingly small scattering angles where the valence- 



quark density is high 



12], nonthermal processes would be more important 



than in the midrapidity region that we are investigating here. Thermal gluons 
will also dissociate the x&(lP) and Xb(2P) states which partially feed the 
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T(1S) ground state in elementary collisions 

Due to the small velocity i; C c of the quarks in the bound state, 
the proper equation of motion for single-particle quarkonium states is the 
Schrodinger equation, with the color-singlet QQ quarkonium potential Vqq. 
Reasonable parametrizations of the potential exist that have been tested in 
detailed calculations of the excited states. 

In particular, the Cornell potential [3] has string and Coulomb part 
Vqq = ar — a e ff/r, where a ~ 0.192 GeV 2 [15] is the string tension, and 
a e ff = 0.471 an effective Coulomb-like coupling constant that accounts for 
the short-range gluon exchange, respectively. 

Although the string contribution to the potential vanishes for light quarko- 
nia in the QGP above T c , it has to be considered at T > T c for heavy 
quarkonia that remain initially confined and are therefore not in thermal 
equilibrium with the plasma. Hence we maintain the string contribution in 
an approximate solution of the gluodissociation problem. 

The string tension of quarkonium decreases with increasing temperature 



T in t 
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ie quark-gluon medium. The screened potential can be written as 



To r 



(1) 



V(r,T) = ar D [l - e~ r / rD ] - 

with r D (T) the Debye radius, r^ 1 = T [47ra s (2iV c -|- iVf)/6p 2 . The number of 
colors is N c = 3, the number of flavors in the QGP taken as Nf = 3, and the 
strong-coupling constant at the T(1S) mass a s ~ 0.2. Because of the inverse 
proportionality of the minimum screening radius that permits a bound state 
to the heavy-quark mass Q, it is much more difficult to dissolve the T(IS') 
in the quark-gluon plasma through screening than the J /if) {IS). 



Table 1: Thermally averaged cross sections < <Jdiss{nS) > in mb for the gluodissociation 
of the T(15'), (2S), (3S) states at four different temperatures T and m g = in 2.76 TeV 
PbPb. The values include screening as described in the text; 2S and 3S states are screened 
completely at high T. 



T 


< <T diss (lS) > 


< a diss {2S) > 


< a diss (3S) > 


(MeV) 


(mb) 


(mb) 


(mb) 


400 


0.094 






300 


0.141 


0.041 




200 


0.124 


0.465 


0.152 


170 


0.080 


0.783 


0.604 



We have calculated the wave functions of the IS — 3S states, as shown in 
Fig. [1] for T = and 200 MeV. They are almost independent on temperature 
for the ground state. For the 2S state, there is an increase of the rms radius 
from 0.50 to 0.77 fm, whereas for the 35* state the rms radius increases from 
0.73 to 1.99 fm. 

Due to the high temperature and ensuing large thermal gluon density 
reached at LHC energies in the midrapidity region, the most important pro- 
cess next to screening that leads to a suppression of Upsilons at LHC is 
gluodissociation. Hence we calculate the gluodissociation cross sections for 
the IS — 35, Xb{\P~)i an d X&(2-P) states as functions of the initial impact- 
parameter dependent temperature in the quark-gluon plasma. Our calcu- 
lation is complementary to the solution of a Schrodinger equation with an 



imaginary-valued contribution to the potential 
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damping of the exchanged gluon as performed in 21J for the Y(IS') and the 
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20| due to Landau 
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Figure 1: (color online) Radial wave functions of the T(15), (25), (35) states (solid, dotted, 
dashed curves, respectively) calculated in the screened Cornell potential for temperatures 
T = MeV (bottom) and 200 MeV (top) with effective coupling constant a e ff = 0.471, 
and string tension a = 0.192 GeV 2 . The rms radii < r 2 y 1 / 2 of the 25 and, in particular, 
35 state are strongly dependent on temperature T, whereas the ground state remains 
nearly unchanged. 
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Xfe(l-P) states. 

The leading-order dissociation cross section of the QQ states through 
dipole interactions with hard gluons (El absorption of a single gluon) had 



221 ] . In an operator product ex- 



been derived by Bhanot and Peskin (BP) 
pansion, they calculate the gluodissociation cross section adi SS with pure 
Coulomb-like momentum eigenstates. This expansion is valid for sufficiently 
small bound-state radii. For an initial gluon of energy E (momentum p) 
the cross section is obtained from the Born amplitude As using the optical 
theorem Im[AB(t = 0)] = E a diss- 

Modifying the BP approach to approximately account for the confining 
string contribution, we use the singlet wave functions computed with Eq.(l). 
Inserting a complete set of eigenstates of the adjoint (octet) Hamiltonian 
—A/rrib + a e ff /(8r) with eigenvalues k 2 /nib (m& ~ 4.75 GeV [15] the bottom 
quark mass) to calculate the dissociation cross sections of the T(1S, 2S, 3S) 
and the x&(lP, 2P) states [23], we obtain 

00 

at s (E) = J dk5 (^ b + e « - E ) l^)' 2 ( 2 ) 



with the wave function overlap integral 

;>oo 

w nS (k)= / drrg s n0 (r)gti(r) (3) 
Jo 

for the singlet radial wave functions 5^ ( r ) °f ^e b quark, and the adjoint 
octet wave functions ^(r). The binding energy of the nS state is e n , and 
the 5 function accounts for energy conservation, k 2 /rrib = E — e n . 

For vanishing string tension a — Y and the corresponding values of the 
binding energy e n , a pure Coulomb IS wave function, and a simplification in 



the octet wave function, this expression reduces to the result in [22| . We can, 
however, evaluate it with the full octet wave function to obtain the T(IS') 
dissociation cross section a\s in terms of the BP expression a bp 



alL ttz(1 + z 2 /4)/ qzy 

1 H — — ) exp \lz arctan q) 



exp (7TZ) 



1 



4 J 



(4) 



with z = l/(4g), and q = y/E/ei — 1. The rhs approaches 1 for z — > 0, 
recovering the BP formula. It approaches for z — > oo, and agrees with the 



result obtained independently by Brambi 



approach in the corresponding limit 



24 



la et al. in an effective field theory 



25| . Their work also considers the 



thermal width of heavy quarkonia due to Landau damping, in addition to the 
break up of a color-singlet bound state into a quark- ant iquark pair that is 
investigated here. In contrast to other assumptions, these authors find that 
breakup is the leading term as compared to Landau damping [25]. 

We obtain new results for the 2S and 3S states from eqs. (2), (3). We 
also calculate the cross sections for the Xfe(l-P) an d Xb(2P) states [23J. The 
gluodissociation cross sections resulting from eqs. 01]), (3) including the effect 
of screening for finite string tension are shown in Fig. [2] for the 15* and 2S 
states. 

One should be prepared to expect modifications in the cross section values 
of the five states from next-to-leading order (NLO) contributions [26], where 
a gluon appears in the final state in addition to the b and b quarks, and hence, 
the phase space is larger than in leading order (LO). However, in {27] it was 
shown that the quasi-free process that corresponds to NLO is less important 
than LO for temperatures T >270 MeV. 

Whereas the heavy quarkonium is not in thermal equilibrium with the 
QGP, it is reasonable to assume that the medium itself is thermalized due to 
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Figure 2: (color online) Gluodissociation cross sections adiss(nS) in mb (lhs scale) of the 
T(15) and T(2S) states calculated using the screened Cornell potential for temperatures 
T =200 (solid curves) and 250 MeV (dotted curves) as functions of the gluon energy E g , 
The thermal gluon distribution (rhs scale, with ra g = 0; solid for T = 200 MeV, dotted 
for 250 MeV) is used to obtain the thermally averaged cross sections through integrations 
over the gluon momenta. 
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Table 2: Calculated minimum-bias suppression factors Raa{1S) of the T(15) state for 
initial central temperatures To = 500-800 MeV in 2.76 TeV PbPb, corresponding results 
for an effective gluon mass of 1 GeV, and ratio of the yields T(2S + 35)/T(15) for effective 
gluon masses m g — (fourth column), and m g = 1 GeV (last column). Suppression factor 
values for finite m g are lower bounds, see text. 



To(GeV) 


RaaO-S) 


RZ(IS) 


T{2S + 35)/T(15) 




0.8 


0.50 


0.61 


0.45 


0.38 


0.7 


0.53 


0.66 


0.49 


0.40 


0.6 


0.56 


0.69 


0.53 


0.44 


0.5 


0.62 


0.74 


0.58 


0.50 



the short equilibration time of about 0.6 fm/c [8], at least in the transverse di- 
rection. Hence, we integrate the gluodissociation cross sections for the IS, 2S 
and 35 states over the gluon momenta p, weighted with the Bose-Einstein 
distribution function of gluons at temperature T to obtain the average dis- 
sociation cross sections for the nS states 

^ „nS ^_ 9d f°° nS P 2 dp , , 

<ad - s> -2^J ad - s[h) exp[E(p)/T}-l (5) 

with E{p) = {p 2 + m 2 g ) 1 / 2 , the gluon degeneracy ^=16, and the gluon density 
as the integral over the distribution function, n g = g<iT 3 ((3)/iT 2 for m g = 0. 
Values for the thermal gluon density at temperatures 170, 200, 300 and 400 
MeV and m g = are n g = 1.25, 2.03, 6.85 and 16.23 fm -3 , respectively. The 
distribution function is shown in Fig. |2] (rhs scale). 

The on-shell gluon energy (p 2 + m 2 ) 1 / 2 is usually calculated assuming 
vanishing gluon mass m g = 0, but we shall also investigate the effect of a 
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finite effective gluon mass, as has 



3ased on lattice QCD results 



29 



3een suggested in quasi-particle models [28] 



30j, with m g ~ 0.5 — 1 GeV. It is argued in 
3l| that the effective mass of the gluons may initially be of the order of the 
gluon saturation scale, m g ~ Q s , which is about 1 GeV at XBjorken = 0.01. 
Here we consider an effect of a finite gluon mass only on the thermal gluon 
distribution, not explicitly in the cross section. The resulting average cross 
section values for finite m g are thus upper limits since a finite gluon mass 
reduces the relative velocity between the T and the gluon. Results for the 
average gluodissociation cross section in mb for T = 170 — 400 MeV, m g = 
are shown in Table [TJ 

The dissociation widths T(nS) of the nS states are then obtained by 
multiplying the average cross sections with the gluon density, T(nS) = n g - < 
a dfss >> an d similarly for the Xb states. To compare with minimum-bias 
data, it is essential to consider the impact-parameter dependence. We as- 
sume a monotonic relation of the initial temperature Tj(6) = T (l — b 2 /b 2 cr ) 
on impact parameter up to a critical value h cr where Ti(b cr ) = T cr = 170 
MeV, with T > 500 MeV, and no gluodissociation beyond b cr , to obtain the 
temperature-dependent suppression factor R(nS, b) (and analogously for the 
Xb{nP) states) prior to feed-down at impact parameter b as 

R(nS, b) = Q[Ti(b) - T c ] exp [-T(nS, b)r{b)} + 6[T C - T t (b)]. (6) 

Here we have used an interaction time of r max ~ 5—8 fm/c in central collisions 



in accordance with hydrodynamic [32J and transport [33| approaches, with 
a monotonic dependence on impact parameter r{b) = T max (l — b/b max ), and 
b max =14.22 fm in PbPb. The numerical values shown in this Letter for 
minimum-bias collisions are obtained after impact-parameter averaging, for 
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Tmax = 8 fm/c. 

This schematic calculation accounts for the essential features of the cen- 
trality on the widths, although a fully time-dependent approach such as 



performed in 



21 



32 



33 



M 



35[ for the destruction of the J '/\& or T me- 



son may yield slightly modified results. It should be noted, however, that 
the gluodissociation widths and in particular, r(lS') have only a relatively 
weak temperature dependence at high temperatures where gluodissociation 
is relevant, such that at a given impact parameter one should not expect 
substantial modifications from a explicit consideration of the dynamics. 

To determine the initial bottomium population vector of all five states 
considered in the cascade calculation, we consider \b populations estimated 
from the CDF feed-down results Q 0.27±0.07(stat)±0.04(sys) for x 6 (lP) -> 
T(15), 0.1 l±0.04(stat) ±0.01 (sys) for Xb(2P) ->■ T{1S); 50.9 % of the IS 
state is directly produced. 

With decay rates for the nS states from the particle data group - includ- 
ing the effect of different branching ratios into the fi + fi~ detection channel 
for n =1, 2, 3 -, we calculate a decay cascade that matches the final popu- 

rn 

lations measured by CMS for pp at 2.76 TeV [10], and thus provides initial 
populations which we use for the PbPb in-medium calculation at the same 
energy. Following the consideration of screening and gluodissociation of the 
five states, we calculate the radiative feed-down cascade in the medium for 
those states which have survived the strong-interaction processes at a given 
impact parameter b, to obtain the final yields in the presence of the QGP. 

Our results for the suppression of the T(IS') state in PbPb relative to pp 
at 2.76 TeV are shown in Table [2] for several initial QGP temperatures T , 
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m g = 0, and 1 GeV. Since the average cross sections for finite gluon mass are 
upper bounds, the corresponding suppression factor values represent lower 
bounds with respect to the influence of a finite gluon mass. 

For initial central QGP temperatures 0.5 GeV < T < 0.6 GeV, our re- 
sults are consistent with the experimental value currently observed by CMS, 
R aa(I S) = 0.62 ± 0.11(stat)±0.10(sys) in minimum-bias PbPb collisions 
3, 111) for both zero and finite effective gluon mass. The suppression factor 
of the ground state may, however, be further reduced by cold nuclear matter 
effects such as gluon shadowing and nuclear absorption. 

For the excited states we take approximate values for the initial popula- 
tions as obtained by CMS for pp at yfi = 2.76 TeV [l0|, X (2S+3S) /T(1S) | 
0.78. Our results for the corresponding population ratio in PbPb are shown 
in the last two columns of Table [2] with m g = and 1 GeV, respectively. 
Even at very high initial central QGP temperatures To they are larger than, 
but still consistent with the experimental value that is currently observed by 



pp — 



icj. 



CMS, T(2S + 3S)/T(lS) PbPb = 0.24 + 0.13/ - 0.12(stat)±0.02(sys) 

As an example, we obtain for To ~ 800 MeV and finite effective gluon 
mass m g = 1 GeV a ratio of 0.38+0.19/— 0.12. Here the estimated theoretical 
error bars account for the uncertainties in the input data that enter our 
calculation. With the presently available data, it seems not yet possible to 
further narrow down the QGP temperature due to the large error bars. Our 
result leaves room for additional suppression mechanisms of the excited T 
states in PbPb collisions. 

The expected physical effect, namely, rising dissociation with rising tem- 
perature, is born out in our approach through the combination of screening, 
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Figure 3: (color online) Suppression factors for the T(1S) state (top) and T(2S + 
35)/T(15) (bottom) calculated in the present work for 2.76 TeV PbPb minimum-bias 
collisions from screening, gluodissociation and feed-down as functions of the 6 = tem- 
perature parameter To for three values of the effective gluon mass m g . The corresponding 
CMS minimum bias results (solid lines) with current statistical and systematic experimen- 
tal uncertainties (dashed lines) are indicated [lOj. The estimated theoretical error bars 
account for the uncertainties in the input data that enter our calculation. 
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gluodissociation, and feed-down, even though the thermally averaged gluodis- 
sociation cross sections first rise and then fall with increasing temperature for 
the 15 state. This is shown in Fig. [3] for the Y(15) and Y(2S+3S)/1S sup- 
pression factors in minimum-bias collisions for three different effective gluon 
masses together with the CMS results Q, and the corresponding error bars. 

To conclude, we have calculated the gluodissociation and screening of 
Y(15), (2S), (35) and Xb states at LHC energies, plus the subsequent radia- 
tive feed-down via the Xb states. The weakly bound 35 state dissolves due 
to screening already at temperatures T > 200 MeV which are close to the 
critical value. For 25 + 35 relative to the 15 state we find a substantial 
suppression due to screening, gluodissociation and feed-down that is consis- 
tent with the value reported by CMS when the experimental error bars are 
considered, but allows for additional suppression mechanisms of the excited 
states. 

We obtain reasonable results for the suppression of the excited T states 
relative to the ground state in PbPb collsions at LHC energies with an initial 
central QGP temperature of 500 MeV< T < 800 MeV, an effective gluon 
mass of m g ~ — 1 GeV, and a central-collision interaction time of r int ~ 5 — 8 
fm/c. Screening and gluodissociation are relevant suppression mechanisms 
in particular for the higher bottomium states. The consideration of the 
subsequent feed-down cascade via the Xb states turns out to be an essential 
ingredient in calculating the suppression of the excited states relative to the 
ground state. 

Although screening of the strongly bound 15 ground state is negligible, 
we find that its gluodissociation is sizeable due to the strong overlap of the 
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IS gluodissociation cross section with the thermal gluon distribution. Its ob- 
served suppression factor Raa(1S) ~ 0.62 in minimum-bias PbPb collisions 
[t| is mainly due to both direct gluodissociation of the 15* state, and to the 
melting and gluodissociation of the XbO-P) an d x&(2-P) states which partially 
feed the IS state in pp,pp and e + e~ collisions. 

For a detailed comparison, one needs data with better statistics that is ex- 
pected to become available from the 2011 PbPb run at the LHC If it turned 
out to be possible to measure the populations of the 2S and 3S states very 
precisely, one could use this as a fairly accurate thermometer for the initial 
temperature To of the quark-gluon plasma. On the other hand, substantial 
deviations from the experimental values might indicate that further mecha- 
nisms contribute to the suppression. It may, however, also turn out that the 
gluon distribution is not fully thermalized, in particular, in the longitudinal 
direction. 
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